Antiferromagnetic spherical spin-glass model 



Danilo B Liarte and Carlos S O Yokoi 

Instituto de Fi'sica, Universidade de Sao Paulo, Caixa Postal 66318, 05315-970 Sao 
Paulo, SP, Brazil 

E-mail: danilo@if.usp.br 

Abstract. We study the thermodynamic properties and the phase diagrams of a 
multi-spin antiferromagnetic spherical spin-glass model using the replica method. It is 
a two-sublattice version of the ferromagnetic spherical p-spin glass model. We consider 
both the replica-symmetric and the one-step replica-symmetry-breaking solutions, the 
latter being the most general solution for this model. We find paramagnetic, spin- 
glass, antiferromagnetic and mixed or glassy antiferromagnetic phases. The phase 
transitions are always of second order in the thermodynamic sense, but the spin-glass 
order parameter may undergo a discontinuous change. 
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1. Introduction 

Mean-field theories for spin glasses have not been limited to their original aim of 
explaining peculiar behaviours presented by some magnetic alloys, but have been applied 
to a large number of complex systems, ranging from biology and optimisation problems 
to information processing [U Ej. The paradigmatic mean- field spin-glass model is the 
Sherrington-Kirkpatrick (SK) model [3]. By the application of the replica method, it 
was found that the low temperature spin-glass phase is described by a solution with 
an infinite number of stages of replica symmetry breaking (ooRSB) according to a 
hierarchical scheme proposed by Parisi [I]. The difficulty posed by the analysis of this 
solution, analytically as well as numerically, encouraged the investigation of simpler 
models that retain some of the essential aspects of the SK model. 

The generalisation from a spin-glass interaction between pairs of spins to an 
interaction among sets of p > 2 spins has attracted considerable interest in this respect 
[5]. For instance, in the limit p — * oo this model becomes equivalent to the random 
energy model (REM) [6], which can be solved with or without the help of the replica 
method. Moreover, the first stage of replica symmetry breaking (1RSB) is shown to be 
exact for this model. Another instance in which the analysis becomes simpler is the 
spherical version of the model [7J [8]. This model can be exactly solved and exhibits 
a spin-glass phase described by a stable 1RSB solution for any p, even to the lowest 
temperatures. 

Experimental works report evidences of a spin-glass behaviour and of coexistence 
of spin-glass and antiferromagnetic orders, both in diluted antiferromagnetic materials 
(e.g. Fe^Mgx.j.C^) [9], [101 |TTJ and in mixed antiferromagnetic compounds (e.g. 
Fe^Mni.^TiOs) [12], [EH]- As a theoretical approach to these systems, the two-sublattice 
SK model was introduced to describe re-entrant transitions from the antiferromagnetic 
to the spin-glass phase [TH, [15] . An extended version of this model was investigated in 
an attempt to reproduce the experimental results observed in the experimental systems 
mentioned previously [T6]. Needless to say, the solution of the two-sublattice SK model is 
hard to be analysed analytically as well as numerically. As a simpler model which retains 
some essential features of the two-sublattice SK model, a two-sublattice version of the 
REM was proposed recently [17] to explain some experimental results in the disordered 
antiferromagnetic system Feo.5Zno.5F2 [Hj. The foregoing experimental and theoretical 
investigations motivated us to consider a two-sublattice version of the spherical spin- 
glass model with multi-spin interactions. The relative simplicity of the model enables 
us to investigate the phase diagrams of the model for the full range of parameters. 
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2. The model 

Let us consider a set of 2N continuous spins distributed in two sublattices, A and B, 
each consisting of N spins. The model is defined by the Hamiltonian 

j N N 

H = ~ J il-irh-jr S h ■ ■ ■ S ir°h " ' ° Or + E 3i<T 3 - H ^2( S i+ °i) > C 1 ) 
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where if is an applied magnetic field and Jo > is the antiferromagnetic interaction 
between different sublattices. The interactions Ji 1 ...i r j 1 ...j r among the set of r spins on 
different sublattices are independent Gaussian random variables with zero mean and 
variance 

/ j 2 ) = M ( 2) 

where the factor (r!) 2 is a matter of convention while the dependence on N is needed 
to ensure an extensive free energy. The spins, Si in the A sublattice and cij in the B 
sublattice, are real continuous variables ranging from — oo to oo. The partition function 
is given by 

where (3 is the inverse temperature and the delta functions impose spherical constraints 
to ensure the existence of a well-defined limit at low temperatures. 

3. The replica approach 

In the replica method the self-averaged free energy per spin is computed by means of 

- (3f = lim lim In (Z n ) r , (4) 

where (• • -)j denotes the average over the Gaussian random variables Ji x ~-i r ^~- 3r []]■ A 
standard calculation leads to the following expression for the free energy per spin, 
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where I = A, B is the sublattice index, a, (3 = l,2,...,n are replica indices and 

(q_r) Q/3 = q^ 13 with qf a = 1. The saddle-point equations for and m" are given 
by 

E W) ^ m p A -(3H + (3J m a B = 0, (6) 
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where 

W = if ~ ™ a A< (8) 

which are coupled to similar equations obtained by the interchange A «-> B. Here, q°f 
and q°g denote the overlap between replicas, 

?f = ^E^ «?=4E°?°f> ( 9 ) 

i i 

and and are the sublattice magnetisations, 

™s = ^E^ < = ^E°?- ( 10 ) 

% i 

To evaluate the free energy explicitly, it is necessary to impose some structure on 
q/ and mf. In this work we consider the replica-symmetric (RS) and 1RSB Ansatze, 
the latter being the most general solution for this model [TJ. 



3.1. The RS solution 



Usually the RS form of the overlap matrix is appropriate for the description of systems 
when there is only a single equilibrium state. We therefore expect this Ansatz to be 
valid in the regions of high temperatures and high magnetic fields. Assuming for each 
sublattice I = A, B, 



ad 

Qi 



;i - qi) 5 a(3 + qi 



mi 



mi, 



the free energy per spin becomes 
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where m; and qj satisfy the saddle-point equations, 
m A 



Qa 



(3H + pJ m B = 0, 



Qa 



In 2n + pJomAmB 



+ (3 2 J 2 rq r A 7 1 <lB = 0, 
coupled to two similar equations obtained by the interchange A 



B. 



(12) 

(13) 
(14) 



3.2. The 1RSB solution 

At low temperatures the glassy behaviour is signalled by ergodicity breaking, with the 
emergence of many inequivalent pure states that are described by the breaking of replica 
symmetry pp. For the two-sublattice model the 1RSB Ansatz takes the form 

qf = (1 " 9u) 5 a/3 + (qu - qoi) e af3 + q i, < = m T , (15) 
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for I = A,B where, 

a/3 J 1 if a and (3 belong to the same diagonal block, , . 

1 otherwise, 

which results in the following expression for the free energy per spin 

2j{1-Qi l ~<li x \1- q u J 

f3 2 J 2 

- In 2vr + (3J m A m B — [1 - (1 - x) q[ A q[ B - x q r 0A q r 0B ] , (17) 

where the saddle-point equations are given by, 

-f3H + (3J m B = 0, (18) 

P J QlA QlB + (l-* A )(l-a.) " TTT^ " °' (19) 
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with, 



<f I = xq p 0I + (l-x)q p 1I , (21) 



coupled to similar equations given by the interchange A <-> B. Moreover, the dimension 
of the diagonal blocks x also contribute the additional equation 



{mj - q i) (q i - qu) , J_ ( 1 - q T \ , 1 f q i - qu 



1 t: 

2 

Note that we have assumed the same x for the diagonal blocks in both sublattices 
because the assumption x A ^ x B leads to the RS solutions described by (jl4til4|) . 



1-qj) x 2 \l-quj x\l-qj 

+^(qUQob-qWib)-0- (22) 



4. The phase diagrams 



We present the results only for r = 3 since the the general features of the phase diagrams 
do not depend sensitively on r. 



4.1. Results for H = 

The zero- field phase diagram is shown in Figure [TJ We found four phases which meet at 
a multicritical point: the paramagnetic (P) phase, antiferromagnetic (AF) phase, spin 
glass (SG) phase and mixed (MX) phase (or glassy antiferromagnetic phase). The P and 
AF phases show no replica symmetry breaking and are described by the RS solution. 
The P phase is characterised by the order parameters 



m A = m B = 0, q A = q B = 0, 



(23) 
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Figure 1. The zero-field phase diagram for r = 3 showing the paramagnetic (P), 
antiferromagnetic (AF), spin-glass (SG) and mixed (MX) or glassy antiferromagnetic 
phases. Solid lines indicate continuous transitions and the dashed lines D1RSB 
transitions. 



and the AF phase by 

m A = -m B ^ 0, q A = qB- (24) 

The SG and MX phases present replica symmetry breaking and are described by 1RSB 
solution. The SG phase is characterised by the order parameters 

m A = m B = 0, q 0A = q 0B < q\ A = qi B , < x < 1, (25) 

and the MX phase by 

m A = -m B ^ 0, q 0A = q 0B < q\ A = qi B , < x < 1. (26) 

If we make the substitutions m A — > m, m B — > —m, q 0A , q 0B — > q and qi A , qi B — > qi 
in (11711221) . these equations become identical to those of one-sublattice p-spin spherical 
model with ferromagnetic interactions —Jo and p = 2r [8J. Thus in the absence 
of an external field the results of the one-sublattice model can be translated to the 
two-sublattice model simply by exchanging the ferromagnetic and antiferromagnetic 
orderings. The P-AF and the SG-MX transitions are continuous and are characterised 
by the appearance of spontaneous sublattice magnetisations m A = —m B ^ in the 
AF and MX phases. The P-SG and the AF-MX transitions are characterised by the 
emergence of 1RSB solution with higher free energy than the RS solution in the SG and 
MX phases. In the SG to the P transition x — > 1 and qu — g / vanishes discontinuously. 
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Thus there is a discontinuous one-step replica-symmetry-breaking (D1RSB) transition. 
The AF-MX transition is also a D1RSB transition to the left of the maximum in the 
boundary of MX phase, but to the right of the maximum g A = Qob and q\A = qiB merge 
continuously at the transition with x < 1. Thus there is a continuous one-step replica- 
symmetry-breaking (C1RSB) transition. We observe that although the spin-glass order 
parameters change discontinuously across the D1RSB transition, thermodynamically it 
is a continuous transition [5l [7] . 
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Figure 2. The field-temperature phase diagram of the model for Jo/ J — 1.5 and 
r = 3. The solid lines represent continuous transitions and the dashed lines D1RSB 
transitions. 



4.2. Results for H > 

The phase diagram of the model in a uniform external field is shown in Figure [2] for 
Jo/ J = 1.5 and r = 3. We again found four phases that meet at a multicritical point. 
The P and AF phases show no replica symmetry breaking and are described by the RS 
solution. The P phase is characterised by the order parameters 

m A = m B ^ 0, q A = q B , (27) 

and the AF phase by 

m A ^ m B, qA ^ qB- (28) 



Antiferromagnetic spherical spin-glass model 8 

The SG and MX phases present replica symmetry breaking and are described by 1RSB 
solution. The SG phase is characterised by the order parameters 

m A = m B ^ 0, q 0A = q 0B , q 1A = q 1B , < x < 1, (29) 

and the MX phase by 

m A ^ m B , q 0A ^ q 0B , q 1A ^ q 1B , < X < 1. (30) 
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Figure 3. The order parameters as functions of the magnetic field for r = 3, 
Jo/ J — 1-5 and T/J — 0.5. The full lines represent m A and wis, the dotted lines 
qaA and q$B, the dashed lines qi A and q±B, and the dash-dotted line represents x. 

The order parameters as functions of the external field is plotted in Figure [3] for 
r = 3, Jo/ J = 1.5 and T/J = 0.5, when the system undergoes transitions from MX to 
SG and from SG to P phases for increasing fields. The MX-SG transition takes place at 
H/ J = 1.32. At this transition the sublattice order parameters m A and m B , qo A and qo B , 
and q\ A and q\ B become identical with x < 1, implying a continuous transition between 
RSB phases. At the SG-P transition at H/J = 3.15 the order parameters q 0A = q 0B 
and q\ A = q\ B merge continuously with x < 1, indicating a C1RSB transition. 

Figure H] shows the order parameters as functions of the temperature for r = 3, 
Jo/ J = 1.5 and H/J = 1, when the system undergoes transitions from MX to AF 
and from AF to P phases for increasing temperatures. The MX-AF transition takes 
place at T/J = 0.58. At this transition the sublattice magnetisations m A and m B meet 
continuously. However q QB and q\ B do not merge continuously with q\ A and q\ B . Also 
x — > 1 at the transition, which indicates a D1RSB transition. The AF-P transition at 
T/J = 0.95 is a conventional continuous transition between RS phases. 
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Figure 4. The order parameters as functions of the temperature for r = 3, Jo/ J = 1.5 
and H/J = 1. The full lines represent tua and itib, the dotted lines qoA and goBj the 
dashed lines q±A and qis, and the dash-dotted line represents x. 




Figure 5. Free energy per spin as a function of temperature for H/J = 1.0, 
Jo/ J — 1-5 and r — 3. The solid and dashed lines represent the 1RSB and RS 
solutions respectively. 



The free energies of the RS and 1RSB solutions for r = 3, Jq/ J — 1.5 and H/ J = 1 
are shown in Figure [5] across the MX-AF transition. We observe that in the MX phase 
the free energy of the 1RSB solution is higher than the RS solution. A similar behaviour 
was also found in the one-sublattice model [7]. In the replica approach this indicates 
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(a) J /J = 0.25 



(b) J / J = 0.5 
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Figure 6. The field-temperature phase diagrams for various values of Jo- The full line 
represent continuous transitions while the dashed and dotted lines represent D1RSB 
transitions. 

Figures (B(a)-(f) show the evolution of the field-temperature phase diagrams for 
increasing values of the antiferromagnetic coupling J . In Figure Eta) for J /J = 0.25 
only the SG and the P phases are present. The P-SG boundary has a maximum. It is 
a D1RSB transition to the left of maximum and C1RSB transition to the right. This 
feature of P-SG transition is common to all the subsequent phase diagrams. In Figure 
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HHb) for J /J = 0.5 the AF and MX phases are present for small magnetic fields. The P- 
AF and MX-SG transitions are always continuous. The MX-AF transition is of D1RSB 
type. In Figure [6]^c) for Jo/ J = 1.0 the MX phase extends to zero temperature. In 
Figure [6]^d) the AF-MX boundary exhibits a C1RSB transition in the low field side. 
This feature of the AF-MX transition which starts at Jo/ J = 1.72 is common to all 
subsequent phase diagrams. In Figure E)^e) for Jo/ J = 3.5 the MX phase is no longer 
present at zero field. Finally, in Figure[6](f) the D1RSB transition line has disappeared in 
the P-SG boundary, and only the C1RSB transition line remains. Also the the D1RSB 
transition in the MX-AF boundary becomes again a C1RSB transition before reaching 
the multicritical point. 

5. Conclusions 

We have studied the thermodynamic properties and phase diagrams of a r-spin 
antiferromagnetic spherical spin-glass model using the replica method. For this class of 
models the first step in the replica symmetry breaking is sufficient [7]. The model is a 
two-sublattice version of the p-spin ferromagnetic spherical spin-glass model [8]. The 
two models become essentially identical in their properties in the absence of an external 
field if the roles of the ferromagnetic and the antiferromagnetic orderings are exchanged. 
The model can also be considered as a spherical version of the antiferromagnetic SK 
model [TU UHl EE] and the antiferromagnetic REM model [17] . 

We have presented a detailed numerical study for the representative case r = 3. The 
phase diagrams comprise the paramagnetic (P) phase, the antiferromagnetic (AF) phase, 
the spin-glass (SG) phase and mixed (MX) or glassy antiferromagnetic phases. All the 
transitions between these phases are continuous in the thermodynamic sense. However 
the spin-glass order parameters may change continuously (C1RSB) or discontinuously 
(D1RSB) in the SG-P and MX-AF transitions. 

Previous studies of the same problem in antiferromagnetic SK model [El [H2, [16] 
and the antiferromagnetic REM model [17] have yielded qualitatively similar phase 
diagrams. However the P-SG and AF-MX transitions are of continuous ooRSB type in 
the antiferromagnetic SK model [Hl[T5j[T6] anc ^ °f continuous C1RSB type in the REM 
model [T7] . 
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